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ABSTRACT 

Hubble Space Telescope observations of the site of the supernova (SN) 2008ax obtained in 2011 and 
2013 reveal that the possible progenitor object detected in pre-explosion images was in fact multiple. 
Eour point sources are resolved in the new, higher-resolution images. We identify one of the sources 
with the fading SN. The other three objects are consistent with single supergiant stars. We conclude 
that their light contaminated the previously identified progenitor candidate. After subtraction of 
these stars, the progenitor appears to be significantly fainter and bluer than previously measured. 
Post-explosion photometry at the SN location indicates that the progenitor object has disappeared. If 
single, the progenitor is compatible with a supergiant star of B to mid-A spectral type, while a Wolf- 
Rayet (WR) star would be too luminous in the ultraviolet to account for the observations. Moreover, 
our hydro dynamical modelling shows the pre-explosion mass was 4-5 Mq and the radius was 30-50 
Rq, which is incompatible with a WR progenitor. We present a possible interacting binary progenitor 
computed with our evolutionary models that reproduces all the observational evidence. A companion 
star as luminous as an O9-B0 main-sequence star may have remained after the explosion. 

Subject headings: supernovae: general - supernovae: individual (SN 2008ax) 


1. INTRODUGTION 

The question of what stellar system produce the differ¬ 
ent types of core-collapse supernovae (SNe) is still open. 
While it is accepted that hydrogen-rich (Type H-Plateau) 
SNe arise fr om relatively low-m ass stars in the red super¬ 
giant phase (jSmartt et al.|2009|) , the nature of hydrogen- 
deficie nt (Types Ib and Ic) SNe progenitors remains un¬ 


clear (Eldridge et al. 2013). One crucial aspect is un- 


derstanding what mechanisms are at play for envelope 
removal. In this respect it is interesting to focus the 
attention on the intermediate group of Type Hb SNe 
(SNe Hb) whose spectra evolve from hydrogen rich to hy¬ 
drogen poor and which has been suggested to originate 
from stars that lost most of their H-rich envelopes. Pro¬ 
posed mass-loss mechanisms are those of strong stellar 
winds or interaction in close binary systems. The former 
mechanism requires a very large mass of the SN progen¬ 
itor, possibly with a Wolf-Rayet (WR) or luminous blue 
variable (LBV) nature. 

Remarkably, the firmest identifications of SN-IIb pro¬ 
genitors in pre-explosion images—i.e., those of SNe 1993J 
and 201 Idh—have been associated to moderate-mass su- 


pergiants of spectral types E to K (Aldering et al. 1994 
Van Dyk et al. 2QQ2[ Mannd et ah |2011| |Van Dyk et 
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al.|2011 ), hotter than the progenitors of Type H-Plateau 
iNe. Moreover, they were suggested to have formed close 

detection of their compan- 


ion stars ( 

Mannd et al. 

2004 

et al.l 2014 

[f. SNe 2UU8ax and 


Eox et al. 2014 Eolatelli 


explosion imaging. The case of SN 2008ax is revised in 
this work using newly available data. Eor SN 2013df , the 


candidate is aga in a cool, K-type supergiant star (Van 


Dyk et al. 2014). These cases appear to indicate a pre¬ 
vailing cTo^se bina ry origin for SNe Hb. However, |Gal- 


Yam et al. (2014) found signatures of a WR-type wind 
in very ealy spectra of the SN Hb 2013cu. 

The SN Hb 2008ax has been an important case of 
study due to its small distance and the availability of pre¬ 
explosion Hubble Space Telescope ( HST) imaging and ex- 
tensive multi-wavelength followup. jGrockett et al. (2008) 
studied the pre-SN images and found a source at tne SIN 
location that was possibly the progenitor. With the then 
available data, the authors selected two out of a series 
of possible scenarios: a WR star that retained a small 
amount of hydrogen, or an interacting binary system. 
The light curves and spectra indicated a relatively low 


20081 IRoming et al. 20091 I'svetkov et al.| 


llVlau- 

rer et al. 20101 Taubenberger et al. 

2011 Jerkstrand et 


indication against a single massive WR progenitor. How¬ 
ever, by inclu ding rotation in the e volutionary models of 


massive stars 


Georgy et al. (2012) found a suitable pro¬ 


genitor born with 20 Mq and ending with the correct core 
mass, hydrogen content, luminosity and color to explain 


the complete set of observations. Groh et al. (2013a) 


reinterpreted the final stage of the rotating model as a 
luminous blue variable (LBV) star and suggested that 
LB Vs may be the progenitors of some core-collapse SNe. 
One intriguing aspect of SNe Hb is the possible re- 
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lation between progenitor radius (or H-envelope mass) 
and circumstellar medium (CSM) density. If the pro¬ 
genitor star can be detected in pre-SN observations, a 
radius can be determined. Alternatively, the shape of 
the UV-optical light curve during the first few days after 
explosion can also indicate the extent of the exploding 
star depending on the stre ngth of the cooling emission 
from the sh o cked envelope ( Woosley et al.||1994 Bersten 


et al. 2Q12[ 2014). The CblVl can be probed through 


radio and X-ray observations and this can be linked to 
the mass-loss history of the progenitor. When a dense 
CSM is present shock-heated gas can also produce opti¬ 
cal emission at late times, typically evident in box-like 
emission lines in the spectra and a slow-down of the light 
curve as compared to the usual radioactive-power decline 
rate. These phenomena were o bserved in th e SNe Ilb 
1993J dM atheson et al.||2QQQa|b|) and 2013df (|Maeda et 
~ 4l[2M^ 

ChevaJ 


al 


Jhevalier & Soderberg 


(2010) suggested that ra- 
di'o properties of core-collapse SlNe provide information 
about the size of the progenitor and proposed that 
SNe lib were divided into those with compact (WR- 
like) and extended (supergiant) progenitors. SN 2008ax 
was placed in the compact c ategory while SN 199 3J be¬ 
longed to the extended one. jBersten et al. (2012), how¬ 
ever, showed that such a connection breaks down for 
SN 201 Idh. Nonetheless, by restricting to the three 
SNe lib with direct progenitor detections and estimates 
of the CSM densit y (SNe 1993J, 2011dh, and 2013df), 
Maeda et al. (2015) suggest that SNe with more extended 
progenitors niay have denser surrounding media, possibly 
due to a larger mass-loss rate near the time of explosion. 
The range of radii goes from ?^200 Rq in SN 201 Idh 
to «600 Rq in SNe 1993J and 2013df. There is a wide 
variety of SNe-IIb light curves observed during the first 
few days after explosion, which may indicate a range of 
progenitor radii. H owever, it is not clea r whether the re¬ 
lation suggested by Maeda et al. (2015) holds in general, 
and if it represents a continuum of progenitor proper¬ 
ties or two distinct groups. Adding more objects to the 
CSM density-progenitor radius diagram is thus of criti¬ 
cal importance to understand the evolutionary paths and 
mass-loss history of massive stars. 

Since the work by Crockett et al. (2008), new HST 


imaging of the field of SN 2008ax was obtained with 
longer exposure times and better spatial resolution than 
the pre-explosion ones. The new images allowed us to 
revisit the SN site in order to study the nature of the 
progenitor. In Section we present the high spatial res¬ 
olution images used in this work and provide estimates 
for the distance and extinction toward SN 2008ax. Sec¬ 
tion is dedicated to the analysis of the pre- and post¬ 
explosion photometry at the SN site. Section presents 
the possible properties of the progenitor object based on 
hydrodynamical and stellar evolution models. Our con¬ 
clusions are given in Section 

2. OBSERVATIONAL DATA 
2.1. High-Resolution Imaging 

This work is based on high-resolution images of the field 
where SN 2008ax went off in NGC 4490. The observa¬ 
tions are summarized i n Table including t he images 
previously analyzed by Crockett et al. (2008) and new 


imaging obtained in 2011 and 2013. The former dataset 
included images in four optical bands obtained at three 
different epochs with the Wide Field Planetary Camera 2 
(WFPC2) mounted on HST. The SN site landed on WF 
detectors for the F450IT, F606IT, and F814IT frames, 
yielding a pixel scale of O'.'l pixel” In the F300IF 
frames the SN site was on the PC chip, with a pixel 
scale of 0'.'046 pixel”The new imaging was obtained 
in two epochs (in July 2011 and October 2013) using the 
Wide Field Camera 3 (WFC3) of HST with a pixel scale 
of 0'.'04 pixel”!. In 2011 the F336VF, F438VF, F606VF, 
F625IT, and F814IT bands were employed, while in 2013 
the imaging was done in the F275IT, F336W, F438IT, 
FbbbW, and F814IT bands. Along with a higher reso¬ 
lution, the new images are deeper than the pre-explosion 
ones. 

Processed HST data were downloaded from the Hub¬ 
ble Legacy Archiv^ (HLA). The WFPC2 observations 
were performed in pairs of exposures with each filter. 
The combined F606IT image had poor quality due to 
a misalignment of both exposures. We thus performed 
the image combination for that band using the DRIZ- 
ZLEPAC package. The WFC3 images from 2011 and 
2013 consisted on four and three integrations per band, 
r espectively. 


Crockett et al. (2008) obtained ground-based, high- 
resolution imaging of the SN while it was still bright in or¬ 
der to accurately determine the location of the SN. These 
observations were carried out in the K band using adap¬ 
tive optics with the Altair/NIRI instrument mounted on 
the Gemini North Telescope. We obtained the raw im¬ 
ages from the Gemini Science Archiv^ and processed 
them through sky subtraction and combination using the 
dedicated tools for NIRI in the gemini package of IRAF. 
The co-added image has excellent quality, with FHWM 
of 0'.'09 (« 4.1 pixels). 

2.2. Distance and Extinction 

Distance and extinction are usually the largest sources 
of uncertainty in the study of SN progenitors from direct 
detections. The distance to NGC 4490 has been com¬ 


puted as 9.6 ± 1.3 Mpc by Pastorello et al. (2008) from 
an average of several measurements. 'I'his value has been 
extensively used in the literature for SN 2008ax. We 
note that the average includes distances derived from 
the recession velocity of the host galaxy. These are 
relatively large distances that involve large uncertain¬ 
ties given the small recession velocity of NGC 4490. 
In this work we instead adopt an updated distance of 
7.77 ± 1.54 Mpc, as provided by the NASA/IPAC Extra- 
galactic Database (NED), b ased on seven measurements 
done with the Tully-Fisher (|Tully||1988[ iTheureau et al 


2007 


Karachentsev et al. and Sosies ( [ibrry et al. 


2002) methods 
Ernst extinction in the direction toward SN 2008ax ap¬ 
pears to be significant. The Milky-Way contribution is 
small, the NED provides E{B — V)mw = 0T19 mag 


based on infrared dust maps of [Schlafiy fc: Finkbeiner 
(|2011|). For the host-galaxy reddening, several estimates 
have been given in the literatu re based on equ i valent 
widths (EW) of interstellar lines. Pastorello et al. (2008) 


® http://hla.stsci.edu 

^ WWW.cade-ceda.hia-iha.nrc-enre.gc.ca/en/gsa 
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TABLE 1 

Summary of HST observations. 


Date 

Instrument / 
Detector 

Filter 

Exposure 

(s) 

Program 

ID 

Program 

PI 

1994-12-03 

Pre-explosion 
WPFC2/WFC F606VF 

160 

GO-5446 

Illingworth 

2001-07-02 

WFPC2/PC 

F300VF 

600 

GO-9124 

Windhorst 

2001-11-13 

WFPC2/WFC 

F450VF 

460 

GO-9042 

Smartt 

2001-11-13 

WFPC2/WFC 

F814VF 

460 

GO-9042 

Smartt 

2011-07-12 

Post-explosion 
WFC3/UVIS F336VE 

1704 

GO-12262 

Maund 

2011-07-12 

WFC3/UVIS 

F438VF 

1684 

GO-12262 

Maund 

2011-07-12 

WFC3/UVIS 

F606VF 

1310 

GO-12262 

Maund 

2011-07-12 

WFC3/UVIS 

F625VF 

1029 

GO-12262 

Maund 

2011-07-12 

WFC3/UVIS 

F814VF 

1939 

GO-12262 

Maund 

2013-10-29 

WFC3/UVIS 

F275W 

2361 

GO-13364 

Galzetti 

2013-10-29 

WFC3/UVIS 

F336W 

1107 

GO-13364 

Galzetti 

2013-10-29 

WFC3/UVIS 

F438VF 

953 

GO-13364 

Galzetti 

2013-10-29 

WFC3/UVIS 

F555VF 

1131 

GO-13364 

Galzetti 

2013-10-29 

WFC3/UVIS 

F8UW 

977 

GO-13364 

Galzetti 


inferred a reddening of E{B — R)host = 0.3 mag from 
the EW of the Na I D l ine an d using the empir i cal re- 
lation ofjTuratto et al.|(|2003). jChornock et ah (2011) 
noted that the IN a I D line is saturated and used the EW 
of the K I A7699 instead to derive E(B — E)hos t = 0- 5 
mag based on the rela tion of |Munari & Zwitter (|1997|). 
Chornock et ah ( 2Q11[ ) also provided an EW of the dif- 
fuse interstellar band (DIB) at 5780 A observed in high- 
resolution spectra. W ith their value of 228 mA and the 
relation presented by Phillips et al. (2013), we obtain 
E{B — E)host = 0.4 ±0.2 mag, assu ming a standard red¬ 
dening law of Cardelli et al. (jl989 ) with Ry = 3.1. 

A more accurate measure of the host-galaxy reddening 
can be obtained from the color evolution. We compared 
the MW-corrected (B—V) colors of SN 2008ax with those 
of a sample of reddening-free stripped-envelope SNe ob¬ 


served by th e Carnegie Supernova Project (CSP; Hamuy 


et al. 


2006). The reddening-free sample is composed of 


SNe with no narrow Na I D absorptions and located 
relatively apart of star-forming regions. An average in¬ 
trinsic {B — V) color curve was constructed from this 
sample between maximum light and 20 days after. The 
intrinsic colors evolve monotonically during this interval 
from {B — V) = 0.35 mag to 1.08 mag, with a disper¬ 
sion of 0.06-0.14 mag. The color curve of SN 2008ax 
reproduced the same shape but shifted toward redder 
colors. We adopt the average shift as the color excess, 
E(5 —E)host = 0.27T0.02 mag. Due to the dispersion in 
the intrinsic colors, we consider a systematic uncertainty 
of 0.1 mag. 

We also compared (V — R) colors evaluated at 10 days 
past V- and R-ban d maximum light with the intrinsic- 
color calibrations ofjDrout et al.|(2011) of ((V — R)}vio = 
0.26 ± 0.06 mag and {{V — R)}mo = 0-29 ± 0.08 mag, 
respectively. Eor SN 2008ax we measure (V — R)ylo = 
0.51 ± 0.03 mag and (V — R)rio = 0.57 ± 0.03 mag after 
correcting for MW reddening. An average of the differ¬ 
ences from the reference values yields E{V — R)uost = 
0.26 ± 0.07 mag. Assuming a standard extinction law 
with Ry = 3.1, we obtain E{B — E)Host = 0-38 ± 0.08 
mag. 

Considering the color excesses derived from observed 
colors, we adopt a host-galaxy reddening of E{B — 


EI hosi, = 0-3 T 0.1 mag . The same value was adopted 
by Crockett et al.| (2008) in their analysis of the progeni¬ 
tor cahHiduTeTTIIthou^^ distance was ~ 20% longer. 

3. ANALYSIS 

3.1. Image Registration and Photometry 

We used the combined NIRI image to find the location 
of the SN on the F8141T images obtained at different 
epochs. We employed 39 stars in common — the same 
stars at all epochs — to align the HST images with re¬ 
spect to the NIRI image. The registration produced x, y 
rms uncertainties of 28,26 mas (0.28,0.26 pixel), 8.8, 9.6 
mas (0.22,0.24 pixel), and 7.2,11.2 mas (0.18,0.28 pixel) 
for the 2001, 2011, and 2013 images, respectively. The 
location of the SN on the F8141T images is shown 
in Eigure and corresponds to pixels [208.07,553.01], 
[512.57,5506], and [3516.25,2246.40] in the WEPC2 
image, the 2011 WCE3 image and the 2013 WEC3 image, 
respectively. We c onfirm the ident i ficatio n of the pre-SN 
object detected by|Crockett et al.|([2008|) in the WEPC2 
images. But the late-time WL'C3 images reveal critical 
information previously unavailable about the possible SN 
progenitor. 

The first thing that strikes the eye from Ei gure[^is that 
the object detected by|Crockett et al.|(|2008|) is iimact re¬ 
solved into four sources in the 2011 irnage. The source to 
the SE matches the location of the SN to within «lcr. We 
identify this source with the fading SN. Indeed, the 2013 
images show that this source has disappeared within the 
noise. We indicate the other three sources with letters A, 


B, and C. The object identified by Crockett et al. (2008) 
as the progenitor candidate is thus contaminated by the 
light from objects A, B, and C. 

Eigurej^ shows the complete set of images used in this 
work. The identification of sources A, B, and C can be 
done in the rest of the WEC3 images from 2011 and 2013, 
except for the bluest bands where the signal-to-noise ra¬ 
tio {SNR) is too low. We performed point-spread func¬ 
tion (PSE) photometry of the pre-SN source present in 
the WEPC2 images, and of the four sources identified in 
the WEC3 images. When the objects were not detected, 
we estimated limiting magnitudes as explained below. 

The pre-explosion WEPC2 images were measured with 
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the DAOPHOT package included in IRAFp^ We em¬ 
ployed the combined frames except in the F606W band 
which showed poor alignment of the individual expo¬ 
sures. In that band we measured PSF photometry on 
the two images and averaged the results. The result¬ 
ing photometry is listed in T able (second row), c om- 
pared with that obtained by jCroSett et al. (2008) on 
the same images (first row). I'he values are consistent 
within the uncertainties. We obtain a larger uncertainty 
in the F606IT band, probably due to a larger inconsis¬ 
tency between the two individual exposures. We thus 
confirm that the magnitude in this band is abo ut 0.3 
ma g brig hter than what was originally reported by |Li et 
al. (2008) who may have used the misaligned image from 
Wel^. 

No object was detected at the SN location in the 
WFPC2/F300IF image (Figure |^. We thus estimated 
a limiting magnitude from the pixel statistics in the re¬ 
gion o f the PC chip around the SN site. Following Harris 
(1990), we defined the limiting magnitude for SNR = 5 
as that of an object with peak pixel value of three times 
the standard deviation of the background (abkg = 0.008 
counts s“^). Due to the lack of point sources in the PC 
chip, we could not compute a suitable PSF that would 
allow us to convert the detection threshold to a magni¬ 
tude. Inst ead, we used the encircl ed energy fractions tab¬ 
ulated by Holtzman et al. (1 995) for apertures of d iffer- 
ent radii, l^rom 'Fable 2(a) in Holtzman et al. (1995), the 
fraction of flux in a radius of 1 pixel versus 120 pixels is 
f 1/120 = 0-33 and 0.38 at the F2551F and F336W bands, 
respectively, with a typical uncertainty of 0.05. We as¬ 
sumed an average was valid for F3001F, and thus the 
limiting flux was Flim(T3001F) = 3crbkg7r // 1/120 = 0-21 
counts s~ ^. Using the flux ca libration of WFPC2 (Ta¬ 
ble 5.1 in|Gonzaga et al.||2010|) the limiting Vega-system 
magnitude corresponding to that flux is miini(T3001F) = 
-2.5 log(Fiim(T300VF) x PHOTFLAM) - Vega^p. With 
PHOTFLAM = 6.137x 10“^^ and Vega^p = 19.406 mag, 
this is miini(T3001F) = 22.8 m ag. This value i s simil ar 
to that of 22.9 mag obtained by Crockett et al. (2008) 

On th e WFC3 image s we used the DULPHOT v2.0 
package ( Dolphin||2000 ). DOLPHOT provides a magni¬ 
tude froni each exposure. We only considered those mea¬ 
surements where the object was detected in all frames. 
In those cases, the final magnitude was obtained as a 
weighted average of the individual measurements. In 
other cases we estimated a limiting magnitude based on 
the local pixel statistics. The detection limit was set 
as three times the background standard deviation, and 
this value was converted into a limiting magnitude in 
the Vega system using the PSF computed for the image. 
The resulting magnitudes and limits are listed in Table 
The photometry of objects A, B, and C remains stal^ 
within the uncertainties between 2011 and 2013, while 
the object identified as the SN in 2011 is no longer de¬ 
tected in 2013. The limits at this latter epoch indicate 
that the SN became, as expected, fainter with time. 

Objects A, B, and C are separated from the SN by 

IRAF, the Image Reduction and Analysis Facility, is dis¬ 
tributed by the National Optical Astronomy Observatory, which is 
operated by the Association of Universities for Research in Astron¬ 
omy (AURA), Inc., under cooperative agreement with the National 
Science Foundation (NSF); see http://iraf.noao.edu 


07105, 07115, and 07113, respectively, as measured in the 
F814W image. At the distance of NGC 4490, these cor¬ 
respond to projected separations of 3.9 ±0.9 pc, 4.3 ±0.9 
pc, and 4.3 ± 0.9 pc. Their profiles are compatible with 
being point sources, although the image resolution with 
full width at half-maximum of FWHM ^ 0708 corre¬ 
sponds to a projected diameter of ~3 pc. Their relative 
brightness and separation from the SN site make them 
blend into one source in the W FPG2 images (see Fig - 
ure fTI . As already indicated by jCrockett et al.| (|2008|) , 
the size of the source in the WFPC2 images is actually 
compatible with it holding more than one star. 

3.2. Spectral Energy Distribution 

In order to remove the contamination of the three neigh¬ 
boring sources from the proposed SN progenitor found in 
the pre-explosion WFPC2 images, we determined their 
spectral energy distributions (SED). For this purpose, 
we transformed the magnitudes into specific fluxes at 
the effective wavelength of each band. The results were 
put in absolute scale by adopting a distance modulus of 
29.42 ± 0.42 mag as given by NED. We further assumed 
the extinction to be the same fo r the four objects as that 
of derived for the SN in Section [2^ The corrected SED 
for objects A, B, and G are shown in the left panel of 
Eigurel^ 

The oED we re fit with stellar atmosphere models from 
Kurucz| (|1993|). The results are shown in Eigure |3| and 
the corresponding effective temperatures and luminosi¬ 
ties are given in Table These fits indicate that the 
objects are compatible with being single stars. The sum 
of the SED from stars A, B, and C was assumed to have 
contaminated the photometry of the pre-explosion ob¬ 
ject detected in the WEPG2 images. The total A±B±C 
SED flux was checked by computing photometry through 
a large aperture of 4.5-pixel radius that encompassed the 
three stars in the 2011 and 2013 images after subtracting 
other nearby sources—including the SN itself in the 2011 
images. The resulting fluxes were comparable (to within 
20%-30%) with the summed SED in all bands. However, 
due to the noise included in the large aperture, the un¬ 
known aperture correction, and the arbitrary aperture 
centering, this aperture photometry was too uncertain 
to be used in the following analysis. We instead com¬ 
puted synthetic photometry using the A±B±C SED in 
the WEPG2 bands and subtracted it from the observed 
SED of the pre-SN object. The resulting uncontaminated 
photometry is shown in the right panel of Eigure and 
listed in Table |2j We note that the contamination ac¬ 
counted for about 34%, 42%, and 55% of the pre-SN flux 
in the F4501F, F6061F, and F8141F bands, respectively. 
The contribution in F3001F was also calculated and it 
was subtracted from the upper flux limit of the pre-SN 
object in that band. 

Figure (right panel) also shows the SED of the re¬ 
maining object at the SN location in the 2011 images 
(at t = 1224 days after explosion), and the detection 
limits from the 2013 images (at t = 2063 days after ex¬ 
plosion). These calculations indicate that most of the 
flux of the pre-existing object—even after removing stars 
A, B and C—has disappeared. The SED at 1224 days 
may be mostly due to the fading SN, as indicated by 
the strong decrease in flux at the F4381F and F814W 
bands observed at 2063 days. The flux limits at 2063 
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TABLE 2 

Revised photometry of the pre-SN object 


Measurement 

FSOOfU 

(mag) 

F450VF 

(mag) 

F606VF 

(mag) 

F814VF 

(mag) 

Reference 

Contaminated^ 

Contaminated^ 

U ncont aminated^ 

> 22.9'' 

> 22.8" 

> 22.9" 

23.66(10) 

23.69(13) 

24.14(22) 

23.36(10) 

23.25(24) 

23.8542) 

22.63(10) 

22.7349) 

23.6I42) 

|Crockett et al.|(|2008|) 
i liis work 

This work 


Note. — Uncertainties given in parentheses in units of 1/lOOth of magnitude. 
^ Photometry of the object detected in the pre-SN WFPC2 images. 

^ Revised photometry after subtraction of contributions from stars A, B and C. 
3cr detection limit. 


TABLE 3 

Photometry from post-explosion images. 


Date 

F275VF 

F336VF 

F438VF 

F555VF 

F606VF 

F625VF 

F814VF 

2011-07-12 


25.836(291) 

SN 2008ax 
26.085(125) 

25.168(056) 

24.877(110) 

24.825(069) 

2013-10-29 

> 25.6 

> 25.7 

> 26.6 

> 26.9 



> 25.6 

2011-07-12 


26.063(326) 

25.529(073) 

Star A 

25.333(123) 

25.257(133) 

24.847(111) 

2013-10-29 

> 25.6 

> 25.7 

25.9844v7) 

25.621(176) 



25.065460) 

Average 

> 25.6 

26.063(326) 

25.613403) 

25.62l476) 

25.333(123) 

25.257(133) 

24.898496) 

2011-07-12 


> 26.2 

25.903(118) 

Star B 

25.223(048) 

25.112(073) 

24.679(061) 

2013-10-29 

> 25.6 

> 25.7 

26.033403) 

25.559(073) 



24.83i403) 

Average 

> 25.6 

> 26.2 

25.954483) 

25.559473) 

25.223(048) 

25.112(073) 

24.719453) 

2011-07-12 


> 26.2 

27.468(396) 

Star C 

25.619(081) 

25.244(135) 

24.348(047) 

2013-10-29 

> 25.6 

> 25.7 

> 26.6 

26.220(131) 



24.379477) 

Average 

> 25.6 

> 26.2 

27.468(396) 

26.220431) 

25.619(081) 

25.244(135) 

24.356440) 


Note. — Uncertainties given in parentheses in milimagnitudes. 


TABLE 4 

Eits to the SED of the putative progenitor 

AND NEIGHBORING OBJECTS. 


Object 

T’eff 

(K) 

L/Lq 

Spectral 

Type 

Star A 

9180 ± 700 

9400 ± 1200 

Al lab 

Star B 

7180 ± 190 

7640 ± 210 

F2 lab 

Star C 

4300 ± 150 

8580 ± 170 

K3 lab 

Progenitor^ 

7600-20000 

2.6-20 X 10^ 

B-A la 


^ One-sigma ranges are given because Tgff and 
L/L© are highly correlated as shown in Figure]^ 


days allow only for a small contribution from any other 
object than the progenitor itself. Such contributions are 
not greater than about 12%, 8% and 18% of the pre-SN 
flux in F4501U, F6061U, and F8141U, respectively. The 
limits in the bluest bands (F2751U and F336W) place 
strong constraints on the near-U V em ission of any puta¬ 
tive left-over object (see Section [tT| . 


4. THE PROGENITOR OE SN 2008AX 
4.1. Progenitor SED 

We fit the revised photometry of the progenitor object 
after removing the flux from the three nearby stars by 
assuming the SED of a s ingle star using the atmosphere 
models of Kurucz (|1^3|). The results of such a fit are 
illustrated in h'igure I'he best-fit model is that with 
Teff = 11200 K and ljOg{L/L q) = 4.69, roughly corre¬ 


sponding to a B6 la star. The uncertainties in the pho¬ 
tometry and lack of detection in the F3001U band allow 
for a relatively wide range of solutions. As shown by 
the confidence region in Eigure temperatures can be 
as low as Teff = 7600 K. Solutions with Teff > 20000 K 
fail to comply with the upper limit in the T3001U band. 
Most importantly, the resulting temperatures and lumi¬ 
nosities are highly correlated. The solutions are within 
the realm of bright supergiants, with spectral types be¬ 
tween B1 and A6. The stellar radii of the fitting models 
lie within a narrow range of 40-70 T©. 

Alternatively, we teste d the atmosphere m odels of SN 
progenitors computed by Groh et al. (2013b) for rotating 
and non-rotating massive stars of solar metallicity. Their 
models predict that SN-IIb progenitors should end their 
evolution as red supergiants, yellow hypergiants, Wolf- 
Rayet or luminous blue variable (LBV) stars, depend¬ 
ing on the initial mass and rotation. In Ei gure right 
panel, we show the SED of the models in Groh et al.| 
(2013b|) computed from synthetic photometry through 
the WFPG2 filters. Based on rotating models and the 
uncorr e cted pre-SN photoni etry from [Grockett et al.j 
(2008), Georgy et al. (2012) showed the progenitor of 
SN 2008ax was well-reproduced by a 20- M© star ini¬ 


tially ro tating at 40% of the critical velocity. Groh et al. 
(2013a) later reinterpreted the final stage of that model 
as an LBV star. At the time such an LBV progenitor was 
able to explain not only the available pre-SN photometry, 
but also it had a relati vely low mass to explain the SN 
light curve (see Section 4.2), and a small amount of H to 
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Fig. 1.— {Top left) Ground-based iF-band image obtained in 2008, while the SN was bright, with Altair/NIRI at Gemini North Telescope. 
{Top right) HST image of the SN site obtained before the explosion with WFPG2 and filter F814VF. {Bottom row) HST images of the SN 
site obtained after the explosion using the WFG3/UVIS and filter F814VF in 2011 {left) and 2013 {right). The four images are shown with 
the same scale of about 3^1 by 2'!Q and aligned with North up and East to the left. The SN location as determined from the registration 
with the Gemini Altair/NIRI image is indicated in each HST frame with a white ellipse of her semi-axis in x and y. The object detected 
in the pre-explosion image is resolved into four objects in 2011, one of which is the progenitor candidate that fades away in 2013. 


produce an SN Ilb. With the revised distance and un¬ 
contaminated photometry presented here, however, this 
model no longer fits the pre-SN photometry, as shown in 
Figure 

A better fit to the revised pre-SN photometry in the 
F450W, FQ06W and F814W ban ds is achieved usin g 
the rotating model with 28 Mq from Groh et al. (2013b). 
However, the UV emission is too bright to coui f^ with 
the upper limit in the F300W band (see Figure An¬ 
other problem of the model is the low H abundance at 
the final stage: the progenitor would be a WN star that 
would produce a Type Ib SN. A further complication is 
the final progenitor mass being ^ 8 Mq, wh ich t oo large 
to explain the light-curve timescale (Section |4.2| ). A pos¬ 
sibly better solution could be that of the non- r otating 
25-Mq model, which according tojGroh et al.| (|2013b ) 
leads to a final configuration that is compatible with a 


SN-IIb event. As shown in the right panel of Figure 
(green squares), the model is slightly underluminous in 
all bands. If, however, the distance is assumed to be 
shorter within the uncertainty the fit to the observed pre- 
SN photometry can be improved. As mentioned before, 
the final mass of the model of 8.2 Mq poses a problem 
when one considers the shape of the SN light curve. 

A different scenario that appears to be recurrent 
among stripped-envelope SNe and SNe lib in particular 
is an interacting binary progenitor. In this case, the pri¬ 
mary (donor) star explodes after transferring most of its 
H-rich envelope to the secondary (accreting) star. As a 
consequence, the secondary may remain after the explo¬ 
sion. The detection limits from 2013 place constraints on 
the brightness of any putative left-over companion star. 
The most luminous companion star allowed is an O9-B0 
type, main sequence star, as shown on the right panel 
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pre-explosion 



WFPC2/F450W 



[li\/FPC2/F814W 



2011 







4/VFC3/F438W 


^WFC3/F606W 


^WFC3/F625W ^ 

A. ^ A 
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/'S 


* ^ >WFC3/F814W 


2013 



Fig. 2.— HST images of the field of SN 2008ax. All images are oriented with North up and East to the left, and the field size is 
approximately F.'G by F.'S. (Left column) Pre-explosion images obtained in 1994/2001 with WFPC2. {Middle column) Post-explosion 
images obtained in 2011 with WFC3/UVIS. {Right column) Post-explosion images obtained in 2013 with WFC3/UV1S. All images have 
been registered to the Altair/NIRI image and the location of the SN is indicated with a 5a error ellipse. The SN is still detected in 2011 
but it is below the detection limit in the shallower images from 2013. Three stars have been identified in the 2011 and 2013 WFC3 images 
that appear blended to the SN progenitor in the coarser-resolution WFPC2 pre-explosion images. Their locations in the images where they 
could be detected are indicated with letters “A”, “B”, and “C”. 
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of Figure In Section [4^ we provide support for a bi¬ 
nary progenitor based on the mass of the pre-SN object 
obtained from hydro dynamical modelling of th e SN light 
curve and expansion velocities. In Section [43| we present 
a possible binary system computed with our stellar evo¬ 
lution code that accounts for the pre-SN properties along 
with complying with constraints from observations of the 
SN itself. 

4.2. Hydrodynamical Models 

A set of explosion models were calculated in order to 
derive the possible physical parameters of SN 2008ax. 
For this purpose we used a o ne-dimensional Lagran gian 
LTE hydrodynamical code (Bersten et al. 2Q11|) ap¬ 


plied to stellar ev olution pre-SN structures (Nomoto 
Hashimoto 1988) with different He core masses. The 


methodology is the same t hat we applied in previous 
analyses of SNe Ilb (e.g., see Bersten et al.[2Qf2 Bufano 


et al.|2Q14D , i.e. we first derived global parameters focus 


ing on the explosion of H-free structures to fit the main 
peak of the bolometric light curve and the expansion ve¬ 
locities, and then we used such parameters to analyze the 
effect on the early light curve of a thin H-rich envelope 
with different extensions. In order to compare models 
and observations on the same time frame, we ad opted an 


explosion date of J D = 2454528.80 ± 0.15 from [Tauben- 
berger et al. (2011). 


I'he bolometric light curve was derived u sing the color- 
based bolometric corrections provided by [Lyman et al.| 
2014). We adopted this approach over integrating the 


ux m the observed bands to obtain a more homogeneous 
result at all epochs independently of the wavelength cov¬ 
erage. The bolometric corrections take into account con¬ 
tributions from the UV and IR ranges that can be sig¬ 
nificant depending on the SN color. We employed the 


BV RI an d g'r'i' light curves published by Pastorello et 
al. (|2QQ8| augment ed by BVRI observations ot jTauben- 


berger et al. | ( |2QiT ). Bolometric corrections were derived 


using the calibrations for (H — R), (H — R), (H — /), 
{V — R), {V — I), {g' — r'), and {g' — i'). The colors 
above were corrected by Milk y-W ay and host-galaxy red¬ 
dening as derived in Section [2^ Observed magnitudes 
were corrected for extinction and the resulting bolomet¬ 
ric magnitudes were transformed to luminositi es u sing 
the distance to NGC 4490 also given in Section [2^ For 
each color a different result was obtained, thus yielding 
seven bolometric light curves. These were finally com¬ 
bined by interpolating each light curve to the epochs of 
the others and averaging the results. 

In order to approximate the photospheric velocities cal¬ 
culated by our models, we adopted line expansion ve¬ 
locities derived from Doppler shifts of the He I A5876 
and Fe II A 5169 absorptions . He I A5876 velocities were 


taken fro m Chornock et al. (2011) and Taubenberger et 
al. (2011), while be II A5169 velocities were obtained from 
the latter work. Roughly after day 10 since explosion 
He I A5876 velocities began to flatten as the photosphere 
receded beyond the bottom of the He-rich layer. At such 
late times Fe II A5169 presumably provides a closer ap¬ 
proximation of the photospheric velocity. At early times, 
the Fe II A5169 velocity has a flat behavior that may be 
due to contamination from other lines as the Fe II lines 
are rather weak. We thus adopted the He I velocities at 
t < 10 days and the Fe II velocities afterward. 


In the first stage of the modeling we derive the ejected 
mass, explosion energy and radioactive nickel production 
by fitting the overall shape of the bolometric light curve 
and expansion velocity evolution. Figure shows the 
observations compared with four models of different pro¬ 
genitor masses: (1) model He3 with He core mass of 3.3 
Mq, explosion energy of 8 x 10^^ erg, and ^^Ni mass of 
0.07 Mq, (2) model He4 with He core mass of 4 Mq, 
explosion energy of 1 x 10^^ erg, and ^^Ni mass of 0.05 
Mq, (3) model He5 with He core mass of 5 Mq, explo¬ 
sion energy of 1.2 x 10^^ erg, and ^^Ni mass of 0.06 Mq, 
and (4) model He6 with He core mass of 6 Mq, explosion 
energy of 2 x 10^^ erg, and ^^Ni mass of 0.07 Mq. 

We consider that models He3, He4, and He5 give a fair 
representation of the observations judging all the uncer¬ 
tainties involved, whereas model He6 provides an overall 
worse fit to the light curve shape and velocity evolution. 
The range of validity of the model parameters is roughly 
given by the three former solutions. Note that mod¬ 
els He4 and He5 provide a better fit on the rising part 
of the bolometric light curve than model He3, but the 
latter gives a better match to the main peak and sub¬ 
sequent decline. While He5 produces a slightly overesti¬ 
mation of the tail luminosity, He4 has a lower peak and 
slightly wider light curve than what was observed. On 
the other hand, He3 gives the best approximation to the 
late-time photospheric velocity evolution as traced by the 
Fe II A5169 line. Even if models He4 and He5 provide 
higher photospheric velocities than the observed ones we 
consider that they are acceptable given the uncertain¬ 
ties involved. Contrarily, model He6 departs from the 
observations in an unsolvable manner. Its light curve is 
too and wide, the initial rise is not well-reproduced, and 
the late-time velocities are too large. If the energy were 
lowered to improve the match to the velocities, the re¬ 
sulting light curve would be even wider and the initial 
luminosity evolution would be even lower. 

This analysis suggests that the progenitor of SN 2008ax 
had a He core mass of < 5 Mq . In order to have this low 
He core mass, the progenitor necessarily had a zero-age 
main sequence (ZAMS) mass below 25 Mq, thus favoring 
a binary origin. Our conclusions are consist ent with pre¬ 
vious results based on analytical modeling ([Pastorello et 


al.|[2008t [Homing et al.[[2009| [Taubenberger et al.[[2UIip. 

Early photometry obtained within a few days from ex¬ 

plosion can serve to reveal the progenitor structure and 
in particular, the ex t ent of its H-rich envelope (e.g., see 


Bersten et al. 2012 2014). During this phase, before 


radioactivity begins to dominate, the light curve is reg¬ 
ulated by the physical conditions of the shock-heated 
ejecta. This in turn leads to an initial luminosity decline 
usually referred to as the “cooling phase”. The more 
extended the progenitor structure, the higher the lumi¬ 
nosity and temperature of the emitting region right after 
shock breakout, and the slower the subsequent decline. 
The effect on temperature makes the cooling phase more 
noticeable in the blue bands than in the red ones. 

Interestingly, SN 2008ax was observed in the UV range 
by the Swif t satellite soon after the explosion (Homing 
et al. 2009). The Swift/WOT data in the UV bands 
showed an initial decline lasting for at least four days 
after explosion. This decline was not observed in the 
optical range, although a slower initial rise may be seen 
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Fig. 3.— Photometry and SED fits of objects in the SN site. (L eft vanet) Stars in the vicinity of SN 2008ax (A: blue squares; B: red 
diamonds; green triangles) with fitted stellar atmosphere SED from |Kurucz| (|1993| (A: solid blue; B: dashed red; C: dot-dashed green). 
The arrows show upper limits in the F275W and F336W bands. Tne black line is the sum of the fitted SED of the three stars. (Right 
panel) Pre-explosion photometry from WFPC2 including stars A, B and C (red squares and red arrow for the upper limit in F300tF), 
total SED of stars A, B and C (black dashed line), uncontaminated pre-SN photometry (blue circles and arrow), SN photometry from 
WFC3 obtained in 2011 at 1224 days after explosion (green diamonds), upper limits of SN flux from WFC3 obtained in 2013 at 2063 days 
after explosion (magenta arrows). The blue dashed-dotted line shows the most-luminous allowed companion star in the case of a binary 
progenitor, corresponding to a O9-B0 type main-sequence star. 




Fig. 4.— Single-star progenitor fits. (Left panel) Confidence region with 68% probability (yellow shaded area) in the HRD from fits of 
Kurucz atmosphere models to the pre-SN object photometry. Red symbols show the parameters of model SED shown in the right panel, 
as labeled. The gray shaded region is forbidden by the detection limit in the F300VF band. (Right panel) Uncontaminated photometry 
of the pre-SN object (blue circles and arrow) compared with model atmosphere SED (red lines). The best-fit model is shown with a solid 
line, rougly corresponding to a B6 la star. The long-dashed line shows the lowest-temperature model within the 68% probability region, 
corresponding to an A6 la star. The highest-temperature model allowed by the detection limit in ^ 30014^ is shown with a short-dashed line 
and represents a B1 la star. Synthetic photometry from model atmo sphe res by|Groh et al.| (|2013b|) are shown with gray squares connected 
by dotted lines. Three of these models that are described in Section [4.1 [ are highlighted: a 28-Mq rotating star (black triangles), a 25-Mq 
non-rotating star (green squares), and a 20-Mq rotating star (red circles). 
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in the Swift b band dRoming; et al. 2009) and in the g' 
band (Pastorello et al.| 2008p. I'hese observations are 
reminiscent ot a post shock cooling phase. We thus in¬ 
vestigated whether the early light curves provided infor¬ 
mation about the progenitor extent based on our hydro- 
dynamical modelling. 

We considered our H-free He5 initial stellar struc¬ 
ture, which has a radius of R = 2 i?Q, and modified 
it by adding thin H-rich envelopes of different extent to 
R = 30, 50, and 100 Rq and masses of 0.02-0.04 Mq. 
With these initial structures we computed the hydrody¬ 
namics using the same explosion parameters as in the 
model shown in Figure The left panel of Figure 
shows the comparison oi the resulting bolometric light 
curves with the observations at early times. Clearly, as 
the radius increases the initial decline becomes more evi¬ 
dent. After about five days since explosion, all the mod¬ 
els converge as the light curve becomes regulated by ^^Ni 
radioactivity. The comparison suggests that progenitor 
radii up to about 30 Rq are allowed. 

A caveat should be made about the accuracy of the 
bolometric luminosity at the earliest times when a sig¬ 
nificant fraction of the fiux is expected to be emitted 
in the UV range. According to Pritchard et al. (2014), 
the observed colors of SN 2008ax indicate that the Irac- 
tion of the total flux included in the Swift UV range is 
well below 10% at all epochs. And even if the bolomet¬ 
ric corrections from Lyman et al. (2014) can account for 
this contribution (see their Section 5.1), any extra flux at 
shorter wavelengths may have been missed. This is why 
in the following we decided to directly compare our mod¬ 
els with light curves in individual Swift/WOT bands. 

For that purpose we assumed a black body emission to 
convert the modeled bolometric light curves into broad¬ 
band light curves. This is a simple assumption, although 
at early times it may provide a valid approximation as the 
ejecta are sufficiently dense. We constrained our analy¬ 
sis to the rt-band data since it showed the most clear 
initial decline. Bluer UV bands showed only upper lim¬ 
its at this early stage. The right panel of Figure shows 
the comparison of our extended He5 models with the u- 
band light curve. Absolute magnitudes were obtained 
from the observations by assuming a d istance modulus 
and extinction as derived in Section [2^ (i.e., i^-band ex¬ 
tinction of Au = 1.47 mag). Although the models fail to 
exactly reproduce the low luminosity at t < 5 days, pre¬ 
sumably due to the black-body approximation, we can 
qualitatively see that the more extended the progenitor 
the larger the initial decline and the later the light-curve 
minimum occurs. Also the extended models produce a 
smaller contrast between the minimum and the subse¬ 
quent maximum of the light curve. The apparent dura¬ 
tion and depth of the observed minimum would indicate 
that the most favored models are those between 30 and 
50 Rq. 

We performed a similar analysis based on the He3 and 
He4 models and found a worse comparison with the ob¬ 
servations. In particular, as the core mass decreases, the 
14-band light curve becomes more luminous. This is a 
consequence of the increased temperatures reached by 
the ejecta due to the shock wave propagation. Although 
the approximations involved may prevent a definitive 
conclusion to be driven, the overall trends of the i4-band 
analysis seem to indicate that the mass of the progeni¬ 


tor was between 4 and 5 Mq and its radius was between 
30 and 50 Rq. Indeed, the location of the pre-SN ob¬ 
ject in t he Hertzsprung-Russell diagram (HRD; see Sec¬ 
tion 4.3) indicates an extended structure for the progen¬ 
itor of R ^ 40 Rq . 

4.3. A Binary Progenitor 

We now consider whether the results from the previous 
sections can be interpreted in the frame of the theory of 
binary stellar evolution. For this purpose we computed 
several models with the aim of accounting for the pre¬ 
supernova photometry and, simultaneously, the observed 
light curve of SN 2008ax. A successful model should pro¬ 
duce a progenitor with a final mass of ^ 4Mq to repro¬ 
duce the SN light curve, and a small amount of hydrogen 
in the outermost layers to satisfy the SN classification as 
Type Ilb. 

We empl oyed the code described in |Benvenuto fc De 
Vito (2003), updated to consider the case of massive bi- 


naries as presented in Benvenuto et al. (2013). Our stel¬ 
lar code handles the mass transfer process m a fully im¬ 
plicit way. We considered models with solar chemical 
composition and ignored the effects of rotation. Further¬ 
more, we assumed that the companion star is able to 
retain a fraction p of the material transferred by the pri¬ 
mary (the one that explodes). P is kept fixed throughout 
the entire evolution of the pair. We ended the calcula¬ 
tions at oxygen exhaustion. Since the remainder of the 
primary evolution occurs during a very short time scale 
and only affecting its core, it is safe to assume that the 
stars do not significantly move any further in the HRD. 

We explored the parameters of the problem—the 
masses of the components and the initial orbital pe¬ 
riod, PUk— based on our study of SN 20IIdh (Benvenuto 
et al.||20I3|). In comparison, for SN 2008ax we found 


slightly larger initial masses and a much shorter 
The observations of SN 2008ax are fulfilled with a pair 
of I8MQ-hI2MQ on an orbit with = 5 days, while 
for SN 20IIdh we proposed a pair with I6MQ-hI0MQ 
and P/^^ = 125 days. 

In Figure we show the resulting evolutionary track 
for the case of {3 = 0.50. Shortly after the end of the 
core hydrogen burning stage, the primary star undergoes 
Roche lobe overflow (RLOF), i.e., a Class B mass transfer 
episode. The main mass transfer episode starts at an 
age of 9.785 Myr and lasts for only 26 Kyr. The star 
detaches from its Roche lobe with a mass of 4.376 Mq; 
thus, the mean mass transfer rate is of 5.08 x 10-^ Mq 
yr“^ while the maximum value is of 2.30 x 10“^ Mq 
yr“^. Core helium burning starts soon after the end of 
the RLOF, and most of it is burned during the large 
blue loop in the HRD. When central helium is exhausted 
the star begins to swell again, evolving to lower effective 
temperature conditions. Then, the primary star begins 
to burn carbon and quickly undergoes the second RLOF 
that lasts up to explosion. The companion star reacts 
to accretion by swelling appreciably, but it does not fill 
its own Roche lobe at any moment; therefore, the binary 
never reaches a contact configuration and thus avoids 
suffering a common-envelope episode. 

The primary explodes at an age of 11.02 Myr. At these 
conditions, the pre-SN has a mass of 4.14 Mq and a ra¬ 
dius of 40.7 Rq. Most of its radius is occupied by the 
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Fig. 5.— (Left panel) Observed bolometric light curve of SN 2008ax (circles) compared with the results of the light curve calculations 
for models He3 (dotted line), He4 (solid line) and He5 (dashed line). (Right panel) Evolution of the photospheric velocity for models He3 
(dashed line), He4 (solid line) and He5 (dash-dotted line) compared with measured line velocities from the He I A5876 (triangles) and 
Fe II A5169 (circles) lines. 




Fig. 6 .— (Left panel) Observed bolometric light curve of SN 2008ax (circles) at early times compared with hydrodvnamical models for 
progenitors with 5 Mq helium cores and different extensions of the H-rich envelopes (lines). Model He5 from Figure [mis that with 2 Rq. 
(Right panel) Observed Swift/UYOT n-band light curve (circles; |Roming et al.|2009|) compared with the same models as in the left panel. 
A black body was assumed to calculate the n-band luminosity from the models. 


outer hydrogen rich layer (outward of 1.4 Rq) containing 
only 0.06 Mq. The total amount hydrogen in the outer 
layer is of ~ 3.6 x 10“^ Mq. At that moment, the com¬ 
panion has a mass of 18.44 Mq with a luminosity and 
effective temperature very similar to those corresponding 
to a 20 Mostar on the Zero Age Main Sequence (ZAMS, 
see Figure^; so, it is somewhat over luminous. 

The parameters chosen for the binary system satisfy 
the effective temperatures and luminosities for both stars 
just before the explosion, and also a primary structure 
capable of accounting the light curve and classification of 
SN 2008ax. However, we should remark that this is not a 
unique solution because the parameters are degenerate. 


For example, we may have chosen a higher initial mass 
for the companion and a lower value of p to reach a sim¬ 
ilar final mass. In any case, it seems virtually impossible 
to find a final mass of the primary of 4-5 Mq —suitable 
to account for the observed light curve—if the initial pri¬ 
mary mass is modified by more than about ±2 Mq. 

The results presented in this section allow us to state 
that binary stellar evolution is able to provide a success¬ 
ful progenitor scenario to explain the main characteristics 
of SN 2008ax. 




















12 


Folatelli et al. 



Log(T.„/K) 


Fig. 7. — Evolutionary tracks in the Hertzsprung-Russell diagram for a binary system of 18Mq+12Mq on an orbit with ~ ^ days. 
The ZAMS locus is shown for reference by a dashed line with some mass values (in solar units) denoted with filled dots. Each track starts 
at the corresponding location on the ZAMS. Thick portions of the primary track indicate stages of active core burning. The conditions at 
the moment of explosion are shown with stars. The primary star undergoes RLOF after hydrogen core exhaustion. When it detaches it 
has a mass of 4.376 Mq and starts a blueward loop on which it burns out its helium core. Since then it swells and undergoes a final RLOF 
on which it completes carbon and subsequent nuclear cycles (neon, oxygen, and silicon burning), and it finally explodes. Simultaneously, 
the companion star accretes half (because we set /5 = 0.50) of the transferred material, reaching a final mass of 18.442 Mq. 


5. CONCLUSIONS 

SN 2008ax has been one of the best-studied SNe Ilb to 
date. Along with extensive, multi-wavelength follow-up 
observations, the existence of deep pre-explosion HST 
images placed this object among a selected group of 
SNe lib for whi ch a progenitor could be characterized, 
as first done by Crockett et al. (2008). Since then, ad¬ 
ditional observations of the SJN site were obtained with 
even higher resolution using the WFC3 camera on board 
HST. The new images allowed us to revisit the pro¬ 
genitor candidate and to discover that it was in fact a 
multiple source. After subtracting the light from three 
nearby stars, the revised progenitor photometry resulted 
significantly fainter and bluer than what was measured 
on the pre-explosion images. Moreover, the new images 
showed a fading source coincident with the SN location 
and demonstrated that most of the light from the pre- 
SN source had disappeared. This observation unambigu¬ 
ously showed the correct identification of the exploding 
star. SN 2008ax is thus the third SN lib, after SNe 1993J 
and 201 Idh, with a firm progenitor detection. 

Comparing the revised photometry with stellar atmo¬ 
sphere models, we found that the progenitor, if single, 
was compatible with a B to mid-A type supergiant star. 
Contraril y, pre-SN ni o dels of WR and LBV stars calcu¬ 
lated by Groh et al. (2013b), some of which had been 
proposed lor SIN 2008ax based on the contaminated pho¬ 
tometry, are not compatible with the new photometry. 


Detection limits from the 2013 images allow for a rela¬ 
tively faint source to be left at the SN site. We estimated 
that any such remaining object could be as luminous as 
an O9-B0 main-sequence star. 

Using hydro dynamical models to reproduce the bolo- 
metric light curve and photospheric velocity evolution we 
arrived at the conclusion that the exploding object had 
a relatively low mass of 4-5 Mq, similarly to what had 
been concluded based on other techniques. From a com¬ 
parison with the early light curves, especially those in 
the UV range, we additionally concluded that the pro¬ 
genitor must have had a relatively extended structure, 
with a low-density H-rich envelope reaching about 30-40 
Rq. Such a low mass and extended radius disfavor a 
massive progenitor that loses the outer envelope via stel¬ 
lar winds and ends its life as a WR or LBV star. We note 
that even if the favored progenitor appears as a blue su¬ 
pergiant star, resembling those of SN 1987A-like events, 
the stellar structure is substantially different in the case 
of SN 2008ax, with a much lower-mass, extended H-rich 
envelope. The difference in H mass is what explains the 
dissimilar light curve shape and spectroscopic properties 
between the SN lib 2008ax and the SN 1987A-hke ob¬ 
jects. 

A viable alternative is that of an interacting binary 
system. In this case, the progenitor is allowed to be less 
massive and the envelope is lost through mass transfer 
to the companion star. We were able to provide such 
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a binary model using our stellar evolution code. The 
model assumed a pair of 18 M 0 + 12 Mq in a close orbit 
(with initial period of 5 days) that naturally explained 
the explosion of a nearly H-free, ~4 Mq star of R ~ 
40 Rq with a remaining main-sequence companion that 
complied with the latest detection limits. The solution 
is not unique but it demonstrates that the scenario is 
plausible. A definitive confirmation and characterization 
of the proposed binary progenitor would require a deep 
search for the companion star in future observations from 
space. 

The case of SN 2008ax sheds light on the origin of 
SNe Ilb, which has implications for the progenitors of 
core-collapse SNe in general. To date there are three 
SNe lib with identified progenitors, namely SNe 1993J, 
2008ax, 201 Idh, and one candidate for SN 2013df. Re¬ 
markably, in all of these cases the progenitors have been 
suggested to belong to binary systems. For SNe 1993J 
and 201 Idh there is further evidence of the presence of a 
hot companion star ( Mannd et al.J|2004 Fox et aL]|2014 


Folatelli et al.||2014 ). This fact points to a possible com- 


mon origin for SINe lib. 

The history of mass loss for SN-IIb progenitors is a cru¬ 
cial aspect to understand their nature. From radio and 
X-ray observations and by studying the late-time optical 
spectra and light curves, evidence of substantial amounts 
of CSM has been found for SNe 1993J and 2013df. Con- 
trarily, SNe 2008ax and 201 Idh appear to have occurred 
in much “cleaner” environments. Such differences may 
indicate a different mass-loss history. In the interact¬ 
ing binary scenario, this may in turn be linked with the 
strength and timing of the mass-transfer episodes and 
with the ability of the a ccreting star to reta in the ma¬ 
terial. As pointed out by Maeda et al. (2015), the cases 
with denser CSM correspond to more extended progen¬ 
itors, of over 500 Rq. With a progenitor radius of only 
~40 Rq as derived from our models, SN 2008ax would 


agree with the same picture. More objects are neces¬ 
sary to confirm if this is a general rule for SNe lib and 
whether the progenitor radius and CSM density follow a 
continuum of properties or if they group into two distinct 
classes. Finally, the physical reason for such a possible 
distinction is yet to be known and would require detailed 
modeling of the progenitor evolution. 
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